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Abstract. - 

We study the Kondo effect in a quantum dot embedded in a mesoscopic ring taking into 
account intradot spin-flip scattering R. Based on the finite-!/ slave-boson mean-field approach, 
we find that the Kondo peak in the density of states is split into two peaks by this coherent 
spin-flip transition, which is responsible for some interesting features of the Kondo-assisted 
persistent current circulating the ring: (1) strong suppression and crossover to a sine function 
form with increasing R; (2) appearance of a "hump" in the independent behavior for odd 
parity, ii-induced reverse of the persistent current direction is also observed for odd parity. 



Intensive attention has been paid to the Kondo effect in a quantum dot (QD), [1,2] in 
which formation of a bound state between the local spin and conduction electrons induces a 
Kondo peak in the density of states (DOS) at the Fermi level Ep, consequently leading to 
strong enhancement of the linear conductance ultimately to the ideal value 2e 2 jh below the 
bulk Kondo temperature T^. Associated with this energy scale, the characteristic length of 
the bound state wave- function, the so-called Kondo screening cloud defined as £g w hvp/T^ 
(vf is the Fermi velocity), is itself mesoscopic size (of the order of 1 pm). 

The easy control of major parameters of these artificial atoms in a wide range makes it 
possible to study the Kondo effect even when the size of the lead is comparable to the size 
of screening cloud Recently, Affleck and Simon [3] suggested that the setup of a QD 
embedded in a closed mesoscopic Aharonov-Bohm (AB) ring is an ideal tool to detect this 
cloud because their perturbation calculations found that the persistent current (PC) in such 
a geometry is a universal scaling function of £r/L (L is the ring circumference), and crosses 
over from the saw-tooth of perfect ring in the limit £j^/L <C 1 to a sine function with vanishing 
magnitude at ^/L ^> 1. Besides, slave-boson mean- field (SBMF) evaluation disclosed that 
this finite size effect affects drastically the Kondo resonance in a way that depending on the 
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total number of electrons (modelo 4). [4,5] There was also a variational calculation about this 
problem. [6] 

On the other hand, the recent issue has been focused on what happens to the Kondo physics 
if an intrinsic spin- flip scattering is stirred in the QD and/or even when the conduction band 
is allowed itself to be ferromagnetic metal. [9-15] By means of the SBMF approach, Ma et 
al. [13] and Lopez et al. [14] found respectively a splitting of the zero-bias peak of the nonlinear 
differential conductance when the spin-flip scattering amplitude R is of the order of the Kondo 
temperature, indicating the Kondo peak splitting and suppression in the DOS. This has been 
further confirmed by using an exact theoretical method, numerical renormalization group 
technique. [15] 

Actually, the electron in the semiconductor QD suffers inevitably intrinsic relaxations 
(decoherence) due to the spin-orbital interaction [7] or the hyperfine-mediated spin-flip tran- 
sition, [8] which can be modeled by a phenomcnological spin-flip scattering term R as done 
in these previous papers. [13-15] Naturally, it would be interesting to study the modifica- 
tion of the Kondo resonance due to the intradot spin-flip transition at a mesoscopic situation 
Q/L ~ 1 and its effect on the PC. 

In this letter, based on the finite-?/ SBMF approach initially developed by Kotliar and 
Ruckenstein, [16] we perform a detailed analysis for this issue. With the help of four slave 
boson parameters, this SBMF approach not only allows to deal with finite on-site Coulomb 
interaction but also is expected to give good results both for the charge fluctuation and spin 
fluctuation of the Anderson impurity model, which enables us to thoroughly analyze this 
problem (the Kondo-type DOS in QD without attaching to leads and the PC circulating 
the ring) sweeping across the whole regions, i.e., Kondo, mixed- valence, and empty-orbital 
regimes. [5, 13] 

We consider a QD with a single bare level ed having a finite on-site Coulomb repulsion 
U and an intradot spin-flip scattering R is embedded in a mesoscopic AB ring, modeled by 
a one-dimensional tight-binding Hamiltonian with the nearest neighbor hopping t, through 
tunneling t n [r\ = L,R). According to the finite-?/ slave-boson approach, [16] we introduce 
four auxiliary boson operators e, p a , and d, which are associated respectively with the empty, 
singly occupied and doubly occupied electron states at the QD, to describe the above physical 
problem without correlation terms in an enlarged space with constraints. In terms of the 
auxiliary particle representation, the effective Hamiltonian for this system (totally N sites 
including the QD) can be written as: 

N-2 

H = -* E E(sWi- + H - c -) + E e ^ da c da 

j=l cr cr 

+UdU + Rc^c dl + Rc^ - E(^44 c icr + tR^c^ N _ la c da z a + H.c.) 

a 

+A«£>t Pff + e t e + d U-l) + Y / ^\c\ a c da vlVa - ^d), (1) 

cr cr 

where c da (cda) is the creation (annihilation) operators of electrons in the QD, the three 
Lagrange multipliers and impose the constraints, and in the hopping term the 
fcrmion operator c& a is replaced by Cd a z a with a many-body correction factor z a = (1 — 
cftd — PvP a )~^ (e.^P a + p\d)(l — e^e — p\p s )~^ . The phase factor (f> is defined by 27r$/$ , 
in which $ and $o = h/e are the external magnetic flux enclosed by the ring and the flux 
quantum, respectively. The coupling strength T between the QD and the ring is associated 
with the ring DOS p(Ep) and hopping t(Ep) at the Fermi energy as it p{E-p)\t{E-p)\ 2 . Here 
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we would like to point out that the intrinsic spin-flip processes lift the level degeneracy of 
the QD, yielding ± R. Because of no special spin direction, it seems that our model is 
equivalent to the one of the QD under external magnetic field (Zeeman effect). However, we 
still prefer the present form since it provides a direct description of spin relaxation in the QD 
and can give us a heuristic comprehension of the spin-relaxation effects on the PC, although 
the spin is flipped without any type of dissipation in this simple model. Moreover, our results 
obtained in the following are also applicable for the case of applying external magnetic field. 
We consider the half-filled case (N e = N) in this letter, and in the continuum limit we have 
E F = 0, p(0) = N/(2nt), \t(0)\ 2 = 2{t\ + t 2 R )/N, and r„ = t%jt (r, = L, R). 

At zero temperature, we utilize the mean-field approximation in which all the boson op- 
erators are replaced by real c numbers (variational parameters). Consequently we can solve 
the noninteracting Hamiltonian Eq. Q exactly by numerical diagonalization and evaluate 
the expectation value of any operators in the ground state |0) by summing over consecutive 
eigenlevels up to the highest occupied single particle energy level (HOEL). Minimizing the 
ground state energy E g [expectation value of the Hamiltonian Eq. Q m the ground state] 
with respect to these variational parameters together with the three constraints form the ba- 
sic self-consistent equations to determine these unknown parameters within the mean-field 
scheme. [5] Once solving the set of nonlinear equations, one can calculate the electron DOS 
p a {oj) for the QD by 

PaH = --\z a \ 2 lm{v\c da {oj + ? :0+ - H)~ l c^\v), (2) 
(\v) stands for the vaccum state) and the PC I by 

1 = "fil* = I E(°l*^ e# 4- lCT ^ - H.c.|0). (3) 

a 

The mesoscopic "Kondo correlation energy" Tk is defined as the binding energy between the 
QD and the ring: 

T K = e d - e F - {E g - E° g ), (4) 

in which eF is the energy of the HOEL corresponding to the Fermi energy and E g is the ground 
state energy of the N — 1 sites tight-binding Hamiltonian with the same electron number N e . 

In the following, we present our calculations only for the mesoscopic ring with £ K /L ~ 1 
(L = N) in the particle-hole symmetric region. We will also take t — 1, = £r = 0.35, F = 
Tl + Tr = 0.245 and choose the on-site interaction U — 8T, for the symmetric (asymmetric) 
case e d = -U/2 (-17/4) corresponding to T£ = 0.013 (0.034) or = 148 (59). 

Figure 1 shows the dot DOS p<t{lo) in the symmetric case ea = —U/2 for different values of 
the intradot spin- flip scattering R in terms of the Kondo temperature X£J = U \f]3 ^exp(— 7r / ' (3) /2tt 
[(3 = —2UY/ed{U + ed)]- It is obvious that in the mesoscopic case L ~ the bulk single 
Kondo peak becomes a set of subpeaks depending on the total numbers of lattice sites N 
(modulo 4) , which are roughly attributed to the patterns of the HOEL and lowest unoccupied 
energy level (LUEL) of the Hamiltonian Eq. (JJJ. [4] For N = An, the energies of the HOEL 
(cf) and the LUEL are both nearly equal to zero, so there is a single main peak in DOS pinned 
at u> — together with several small peaks located at other single particle energy levels. For 
N = An + 2, ep is negative with a distinguishable value and LUEL is equal to — ep. As a 
result, two Kondo peaks appear around lo = with the same amplitudes. While in the cases 
of odd parity (N — An ± 1), ep is nearly below or above zero, respectively, leading to a single 
Kondo resonance slightly shifted from the Fermi energy. This shifting of the Kondo peak from 
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Fig. 1 - The dot DOS for the symmetric case td = —17/2 at various spin-flip processes R without 
magnetic flux. The number of the lattice sites N are (a) An, (b) An + 1, (c) An + 2, and (d) An + 3. 
Inset: p CT (w) around uj = at TV = 4n. 

Fig. 2 - The dot DOS for the asymmetric case = —{7/4. The other parameters are the same 
as in Fig. 1. 



the Fermi energy ep in the odd parity will cause the reduction of PC magnitude (see Fig. 4 
below). On the other hand, introduction of the spin- flip scattering R always reduces and splits 
every Kondo peak regardless of parity. Moreover, with increasing R the DOS suppresses in 
a more pronounced way, and the two split Kondo peaks are entirely symmetric for N = An. 
However, the situation is somewhat different for the asymmetric system = — {7/4 as shown 
in Fig. 2. The profiles of DOS wholly move a little bit towards the high energy. This is 
because the asymmetric system allows charge fluctuations to certain extent and consequently 
increases the resulting single particle energy levels delicately. Associated with the spin-flip 
induced splitting, this feature results in a novel i?-dependent behavior of the PC (see below). 

We also study the splitting width 8 of the main Kondo peak at N — An and plotted it as a 
function of R in Fig. 3(a), which shows a significant reduction of the width 5 in contrast to the 
predicted value 2R by measuring nonlinear differential conductance through a QD connected 
to macroscopic leads. [14] Actually, the exact Bethe-ansatz solution for the QD also reveal 
explicitly this large suppression of the peak splitting due to the strong correlations in the 
presence of the external Zeeman field. [17] Besides, Fig. 3(b) displays that, for the fixed R, the 
peak splitting shrinks when the QD departs away from the particle-hole symmetric point. In 
Fig. 3(c) and Fig. 3(d), we summarize the calculated mesoscopic Kondo temperatures 7k with 
variation of R by using Eq. Q. We find that: (1) The finite size effect remarkably increases 
the binding energy between the QD and the ring, so the mesoscopic Kondo temperatures 7k 
are much higher than the corresponding bulk values 7^; (2) Charge fluctuations can decrease 
this binding energy (because it increases the single particle energy levels), and accordingly, 
lead to a smaller 7k in the particle-hole asymmetric region; (3) The behaviors of 7k with 
rising of R depend on the parity: For even parity N = An ± 2, 7k decreases first with small 
amplitude and then begins to increase, while in the case of odd parity N — 4n±l, 7k increases 
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Fig. 3 - (a) Splitting S of the main Kondo peak vs. R for ea = — U/2 and —(7/4 at N — An and zero 
magnetic flux, (b) 8 as a function of td at R = 0.5F. Kondo temperature Tk/7k vs - R f° r ( c ) the 
symmetric and (d) the asymmetric systems at mesoscopic condition L ~ £k and cj> = w/2. 



monotonously, implying that the Kondo screening cloud shrinks due to the spin-flip scattering. 




Fig. 4 - Normalized PCs vs. magnetic flux at various R for ed = — U/2 (a)-(d) and ed = — U/4 (e)-(h). 
Jo = 2t/N is the corresponding PC of an ideal N sites ring. 



Now we turn to calculation of the PC circulating around the ring, which is believed to be 
a useful way to detect the Kondo resonance due to the fact that the PC is closely related to 
electron transmission through the QD embedded in the ring. We illustrate in Fig. 4 the in- 
dependent PC I/Iq vs. magnetic flux at mesoscopic condition for ed = —U/2 and e d ~—U /4. 
In the absence of spin-flip transition, (1) the PC has an analogous structure of an ideal ring, 
implying high transmission through the QD; (2) its direction depends on the parity of the sites 
number N; (3) the magnitude of the PC for the odd parity is much smaller than that for the 
even parity because of the Kondo peak slightly shifting in the case of the odd parity; and (4) 
the perfect symmetric structures for the particle-hole symmetric point ed — — U/2 guarantee 
I 4n+1 ((j)) = J 4 ™ _1 ((/)), but there are no this feature for the asymmetric systems. As expected, 
applying spin-flip scattering generally lowers the magnitude of the PC due to weakening of 



6 



EUROPHYSICS LETTERS 




Fig. 5 Fig. 6 



Fig. 5 - PC vs. energy level td/U of the QD for different parities and several spin-flip transitions at 

(f> = 7T/4. 

Fig. 6 - The PC as a function of spin-flip scattering R at different energy levels of the QD. 

the Kondo resonance and varies its profile gradually to a rough sinusoidal, even though the 
screening cloud is also reduced at the same time. More strikingly, strong spin-flip scattering 
even changes the sign of dl 4:n+1 (4>) / d(f> at small <p for the asymmetric system from positive 
to negative, indicating that the symmetry of the HOEL changes and the system becomes 
diamagnetic. 

Figure 5 displays the PC as a function of the QD level. For even parities, the PC exhibits 
symmetric feature around the central point e<j = — U/2 and unique magnetic response prop- 
erties (being paramagnetic for N = An, while diamagnetic for TV = An + 2) in the whole range 
of energy levels. In the region where the Kondo effect plays a dominating role, the PC is 
obviously enhanced in contrast to its corresponding ideal value {I / h > 1). However, signifi- 
cant suppression of the Kondo resonance due to increasing R can reduce the Kondo-assisted 
PC rapidly. But in both non-Kondo and deep Kondo regimes, we observe this reduction gets 
gradually weak and even a slight increase because of charge fluctuations. Finally at a large 
scattering R = LOT, the PC shows a Coulomb blockade pattern in analogy with the behavior 
of linear conductance through QD over the Kondo temperature. For the case of odd parity, the 
situation is quite different. The Kondo-assisted PC is asymmetric and a remarkable feature 
can be clearly seen I 4n+1 (ed/U) = I 4n+3 (\ed/U\ — 1) due to the symmetry properties of DOS 
in the finite-l/ Anderson model. Two more striking results can be observed (we only address 
for TV = An + 3 here): (1) The PC is a monotonous function of the spin-flip R at > and 
ed < —U/2, whereas there is no universal R behavior in the middle widow; (2) As mentioned 
above, the PC may change direction with increase of R for < —U/2. We summarize these 
independent features of the PC more clearly in Fig. 6, in which we plot the PC as a function 
of R. Explicitly in Fig. 6(d), the PC exhibits a "hump" structure for e d /U = -\ and -± 
as increasing R, which is due to the fact that the Kondo resonance is shifted away from the 
Fermi energy as shown in Fig. 2(d). 

In short, we have investigated the Kondo resonance in a QD embedded in a mesoscopic 
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AB ring taking into account intradot spin-flip scattering R or application of magnetic field 
in the QD by means of the finite-l/ SBMF approach. The Kondo resonance splits into two 
clear peaks (Zeeman effect) with a largely reduced width when increasing R. We also reported 
the i?-dependent PC behaviors in both Kondo and non-Kondo regimes, which can measure 
this Kondo effect in an equilibrium condition. This model is experimentally accessible for the 
update techniques and provides a possible test for the influence of spin-flip transition on the 
Kondo screening. 
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